Ordered valence bond states in symmetric two-dimensional spin-orbital systems 



Guang-Ming Zhang^'^ and Shun-Qing Shcn^ 
^Department of Physics, the University of Hong Kong, Pokfulam Road, Hong Kong, China. 
^ Center for Advanced Study, Tsinghua University, Beijing, 100084, China 

(February 1, 2008) 

We consider a superexciiange Hamiltonian, H = — X]<ij>(2Si ' ~ |)(2Ti ■ Tj — i), wliicli 
describes systems witii orbital degeneracy and strong electron-piionon coupling in the limit of large 
on-site repulsion. In an SU(4) Schwinger boson representation, a reduced spin-orbital interaction is 
derived exactly, and a mean field theory has been developed by introducing a symmetric valence bond 
pairing order parameter. In one dimension, a spin-orbital liquid state with a finite gap is obtained. 
On a two-dimensional square lattice a novel type of spin-orbital ferromagnetically ordered state 
appears, while spin and orbital are antiferromagnetic. Moreover, an important relation has been 
found, relating the spin and orbital correlation functions to the combined spin-orbital ones. 
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There have been much interests in the properties of 
Mott insulators with orbital degeneracy where the 

electron configuration has an orbital degeneracy in addi- 
tion to the spin degeneracy. Due to the interplay between 
the spin and orbital degrees of freedom, a rich variety 
of spin and orbital ordering effects have been displayed 
p-pl , and a new combined spin-orbital degrees of freedom 
may introduce some new physics to the transitional metal 
oxides [p^-p^. Recently, in order to describe the low- 
energy physics of an insulating crystal with one-electron 
per site with double orbital degeneracy in the limit of the 
large on-site Hubbard repulsion, a superexchange Hamil- 
tonian was proposed | |T^ , p^ 

= - ^ (2S,.S,-i)(2T,.T,-i), (1) 

<ij> 

where two isotropic Heisenberg antiferromagnets (AFM) 
coupled by a quartic term on equivalent bonds, and Si 
and Ti denote the spin-1/2 and orbital-1/2 operators 
at a lattice site i, respectively. The physical condition 
to derive the above effective model Hamiltonian is that 
among the possible two-particle states obtained upon vir- 
tual hopping, the inter-orbital singlet is the lowest energy 
state due to the dynamical Jahn- Teller effect. The con- 
dition might be realized in the large body of new molec- 
ular compounds based on Cgo |l8| or layered fuUerides 
and some two-dimensional co poly mers, such as vinyli- 
dene fluoridetrifluoroethylene |l9|. Based on numerical 
calculations, Santoro et al. argued that the ground state 
of the model shows a spin-Peierls-like dimerization in one 
dimension (ID) while on a two dimensional (2D) 

square lattice no evidence of order of any kinds had been 
found that the ground state is concluded to be a spin 
liquid of resonating valence bonds (VB) [|l^. In order 
to put these results on a firm ground, it is thus of great 
interest to develop new approaches to this model. 

In this Letter, we use an SU(4) Schwinger boson rep- 
resentation 1^ to denote both the spin-1/2 and orbital- 
1/2 operators at the same time, and the model Hamilto- 
nian is simplified to a reduced interaction, describing an 



symmetric pairing attraction among the nearest neigh- 
bor hard-core bosons. By introducing an symmetric VB 
order parameter, a mean field theory is developed. In 
ID, a spin-orbital liquid state with a finite gap in the ex- 
citation spectra is obtained, corresponding to a quantum 
disordered VB state. On a 2D square lattice, due to the 
presence of the Bose-Einstein condensations we found a 
novel spin-orbital ferromagnetically (FM) ordered state, 
corresponding to a short-ranged VB crystal state, while 
both the spin and orbital degrees of freedom form an 
AFM ordering. 

The model Hamiltonian can be rewritten as follows 

H = 2Y, (S. -S, +T, -Tj) 

<'*j> 

- 5] (2S..S, + i)(2T,.T,-|-i), (2) 

<i,]> 

where the direct quadratic couplings among spins and or- 
bitals are antiferromagnetic, respectively, while the quar- 
tic coupling between spins and orbitals is ferromagnetic. 
First of all, the model Hamiltonian has an SU {2)(E)SU{2) 
symmetry, representing rotational invariance in both spin 
and orbital space, and also interchange symmetry be- 
tween the spins and orbitals. Moreover, one also notices 
that the total spin, orbital, and combined staggered spin- 
orbital operators 

3 3 3 

where Q is the AF reciprocal vector, generate an SU(4) 
Lie algebra, and also commutes with the model Hamil- 
tonian. Therefore, Eq.(l) possesses the SU(4) symme- 
try |l^. For a state with the SU(4) symmetry, implying 
that the state is invariant under the rotation in the SU(4) 
space, we have an identity for the static correlation func- 
tions 



1 



where (• • •) represents the expectation value on the SU(4) 
symmetric state. The underlying physics for the iden- 
tities is very clear: if the spin and orbital degrees of 
freedom are antiferromagnetically correlated, the com- 
bined spin-orbital correlation functions will be ferro- 
magnetic. This property originates from the additional 
phase factor in the generators L"^. It should be men- 
tioned that another spin-orbital model Hamiltonian H' ~ 
Y.<i,j> (2Sj • Sj + i)(2T, -Tj + 5), the natural generaliza- 
tion of the SU(2) Heisenberg spin model and exactly solu- 
ble in ID, has also been used to understand some orbital 
properties of transition metal oxides ||l0|-|l5[]. However, 
although these two models have the same symmetry, the 
fifteen generators of the SU(4) symmetry group are not 
the same, and the identities of the correlation functions 
induced from this symmetry are also different. In this 
sense, the physics involved in these two models are in- 
dependent, leading to different spin-orbital properties of 
transitional metal oxides. 

As we know, the Schwinger boson mean filed approach 
has been a successful theory in describing the low-energy 
excitations of the conventional FM and AFM Heisenberg 
spin superexchange models In particular, start- 

ing from short-ranged VB order parameters, it can pro- 
duce either quantum disordered or quantum long-range 
ordered magnetic states. In the present model when two 
sets of SU(2) Schwinger bosons are introduced to de- 
note the spin-1/2 and orbital-1/2 operators separately, 
we have 

2S, • Sj + ^ = ^ 4«"'./3"]/3"i^" 

a,f3 

2Ti • Tj + — = dl^^^di^ndj j^dj^rn 

■m,n 

with local constraints '^la'^^a = 1 ^'^'^ Em = 
1, where the indices a, (3 and to, n are taken values 1 
and 2 corresponding to eigenstates of = ±1/2 and 
= ±1/2, respectively. Both ai^a and di^m operators 
satisiy the boson commutation relations. To treat the 
quartic spin-orbital interaction term on an equal footing 
as the same as the spin-spin and orbital-orbital superex- 
change interactions, the fact has been used that in the 
model Hamiltonian the Hilbert space on each lattice site 
consists of four basic states in terms of \S^\T^): 




These four states form a set of local basis to represent 
the SU(4) symmetry group. The conventional SU(4) 
generators S'^[i) act on a basic state |r/)i according to 
the equation S^{i)\rj)i = 5^^ri\^)i with a local constraint 
Sj^{i) = 1, where the indices ^ and v are taken val- 
ues from 1 to 4 corresponding to the four eigenstates of 



|±^; ±5). The SU(4) Lie algebra is defined by the com- 
mutation relation 

In terms of four-component hard-core bosons, we have 
= with the local constraint b\,^^i,ti — 1- 

Therefore, the quartic spin-orbital exchange interaction 
can be written as 

(2S, • S, + i)(2T, ■T, + \)^Y. f^lAA'^h,' (4) 

By a projection procedure, both the spin-spin and 
orbital-orbital quadratic superexchange interactions can 
also be expressed in terms of the four-component hard- 
core bosons. Fortunately, most of the resulting terms are 
found to be exactly canceled with the quartic spin-orbital 
exchange interaction, and the resulting Hamiltonian is 
reduced to a simple and compact form 

<i,3> 

X [(fej,46j,i + bj^ibi^4) - [bj^zbi,2 + 6j,26i,3)] , (5) 

Actually, this reduced model Hamiltonian has explicitly 
displayed that the main physics of the symmetrically cou- 
pled spin-orbital interaction is to induce a unique attrac- 
tive pairing instability among the nearest neighbor hard- 
core bosons. We would also emphasize here that the quar- 
tic spin-orbital exchange interaction plays the same role 
as that of the spin-spin and orbital-orbital superexchange 
terms, as both quadratic and quartic exchange interac- 
tions have similar forms and are treated on the equal 
footing. 

In the reduced spin-orbital model Hamiltonian, the lo- 
cal attractive interaction shares closely resemblance with 
the effective model Hamiltonian in the conventional BCS 
superconductivity theory. Such a feature allows us to 
introduce a unique short-ranged VB pairing order pa- 
rameter 

^ = -{[bjAbi,i + bj^ibiA) - {bj,3bi,2 + ^^,2^1,3)), 

which can assume to be real. A mean field theory is 
followed to develop naturally, leading to the mean field 
Hamiltonian 

H,nf = A ^ bl^b.f, -XN + ZA'^N 

+^ E [(blAA + blA.i)-(f'lA3+blA2) + ''-^- 

<i,3> 

where the local constraint has been imposed on average 
through a Lagrangian multiplier, N is the total number 
of lattice sites, and Z is the number of the nearest neigh- 
bors. In the momentum space, Hmf becomes 



2 



mf 



k,/j, 



(6) 



where 7k = | 



„ik-(5 



and the vector 5 points to the 



nearest neighbor sites. When a Nambu spinor is defined 

by^''l'(k) = (6^,6[2'^L3'4,4;^-k,4:^-k,3,&-k,2,fc-k,l), 

the mean field Hamiltonian is expressed in a compact ma- 
trix form 



1 



k 



with i?m/(k) = Af^o + 2ZA7kl^2- 
Lagrangian is thus given by 



iXN, (7) 
The corresponding 



Here the generalized 8x8 Dirac matrices have been 
defined by Oi = az ® (To ® ao, O2 = ctx ® Oz ® 
where ax, cry, cr^, and <jq are three Pauli 2x2 ma- 
trices and unity matrix, respectively, fio is the 8x8 
unity matrix, fiiand 1^2 obey the anticommutation rela- 
tion [r^i, ^12]+ = 2ilo. From the Lagrangian, the bosonic 
Matsubara Green function is given by 



G{\i,iuj„ 



'iujn^i — Af^o + 2ZA7kri2 
+ [A2 - (2ZA7k)2] 



(8) 



which implies that the bosonic quasiparticle excitations 
form a continuum band, and their dispersion relation is 
a;k=\/A^ — (2ZA7k)^. Moreover, the free energy can be 
evaluated as 



mf 



-El- 



k,M 



2 sinh( 



/3^^k^ 



ZA^TV - ?,XN, (9) 



and the saddle point equations are derived as follows 
2A 



-E 



-E 

k 



VA2 - (2ZA7k)2 

(2^7k)^ 
VA2 - (2ZA7k)2 



[2r7,B(ci;k) 
[2ns (wk) 



1] =3, 
1] =Z. 



(10) 



The first equation corresponds to the request of the local 
constraint on average J^fii^l fj,biti) = 1, the second one is 
to determine the short-ranged VB order parameter self- 
consistently. 

To get the collective excitations of the model, the dy- 
namical correlation functions of the spin, orbital, and 
combined spin-orbital density operators should be calcu- 
lated. The corresponding density operators can be writ- 
ten in terms of the Nambu spinor introduced above, 



Tt = ivl/t(r^)f^4vl/(r,), 
Lf = i*t(r,)a,vl/(r,), 

where fl^ = (Tz®cfQ®(Jz, — cFz®(Jz®cFi), and — a^® 
<Jz ® Tz are introduced. Except for the anticommutation 
relations [rii,ri2]_|_ = [f^2,f^3]_|_ = [f^2,f^4]+ — 2r2o, all 
other relations between the Dirac matrices from f2i to 
fis satisfy the commutation relations. 

In ID system, we have Z = 2 and — cos fc. At T = 
OK, there are no particle excitations nsi^i^k) = 0. When 
the summations over momenta are converted into the in- 
tegrals within the first Brillouin zone, the self-consistent 
equations are easily solved, leading to A ~ 1.67476 and 
A ~ 7.31391. The corresponding ground state energy per 
site is given by £g ~ —2A^ ~ —5.60964, and there is a fi- 
nite gap in the bosonic quasiparticle excitation spectrum, 
Wo = VA^ — 16 A^ ~ 2.93533 at the momentum k — 
and k — TT. These results have shown that the ground 
state of the model is a quantum disordered state of res- 
onating valence bonds, in qualitative agreement with the 
previous numerical calculations p^ . Moreover, the dy- 
namical susceptibilities at zero temperature can be easily 
evaluated as follows 



dk [A^ ± 16A2cosA: cos(fc + q) 

16 [ UJkl^k+q 

1 1 



LU + LUk + LJk+q + if] 



Uk - LUk+q + ir] 



- 1 
, (11) 



where rj is a positive infinitesimal, both the spin and 
orbital dynamic susceptibilities are equal to each other 
and choose the lower "—"sign, and the combined spin- 
orbital dynamic susceptibility takes the upper sign. 
Then, the dynamic density spectra are easily obtained 

Imxs(g ^ = ImXT(g ->Q,uj)^ ImXL(g ~^ 0,w) 

^ ■ ~ (12) 



\uj\ \ uj^^ 4A2 + 64A2 



where 2\/\'^ - 16A2 < \u!\ < 2A. Clearly, there is a fi- 
nite gap Er = 2cjo — 5.87066 in both the spin and or- 
bital collective excitations at q = Q, and in the combined 
spin-orbital collective excitation at g = 0. Therefore we 
conclude that the ground state is a quantum disordered 
spin-orbital liquid state of short-ranged valence bonds 
with a finite gap in their collective excitations. 

On a 2D square lattice, Z — A and 7k — ^{coskx + 
cos ky). At zero temperature, the conversion from the 
summations over momenta to the integrals will be in- 
valid as (8A) A. Following the analogous treatments 
for the Bose- Einstein condensation ||2^ , we separate the 
divergent terms at k*= and k*= Q from the summa- 
tions to yield 
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(27r)2 ^ 
dPk 



- (8A7k 



4^5 (Wk- ) 



(27r)2 - (8A7k)2 ^1 - (SA/A)^ 



3 

2' 

4^ <«' 



When (8A) A, the boson condensation happens, and 
the ratio p ~ 4nB(tJk*)/\/l ^ (^A/A)^ is finite, which 
can be determined from the first equation by setting 
A = 8A inside the integrah p ~ 0.107. Similarly, from 
the second equation one can also obtain the short-ranged 
VB order parameter A ~ 1.3159 and thus A ~ 10.5271. 
Moreover, the ground state energy per site is also eval- 
uated to be Eg = — 4A^ ~ —6.9263. The bosonic quasi- 
particle excitation spectrum now becomes linearly de- 
pendence of momentum near the minimal points k* = 
and k*= Q. When the boson condensation is carefully 
treated, the dynamical structure factors for spin, orbital, 
and combined spin-orbital collective excitations can also 
be calculated through the fluctuation dissipation theo- 
rem. After some algebra, the final results are given by 



^s(q -> Q, c^) = 5T(q -> Q, ^) = ^L(q -> 0, oj) 
= — ^Ms.gn.-f-K 




1, (14) 



where K(a;) is the complete elliptical function of the sec- 
ond kind. In the range of < \io\ < 2A, there is a 
sharp resonance in the collective excitations at zero fre- 
quency, corresponding to the formation of a novel long- 
range ordered spin-orbital crystal state of resonating va- 
lence bonds. In such a state, both the spin and orbital 
degrees of freedom form an AFM long-range order sep- 
arately, while the combined spin-orbital degrees of free- 
dom has an FM long-range ordering. The corresponding 
magnetizations can be read off from the coefficient of 
the delta function m — p/y/Sn ~ 0.021, which is much 
smaller than that of the SU(2) antiferromagnetic Heisen- 
berg spin model. Actually, this result is consistent with 
the identity in Eq.(|3|), as the VB state in the mean field 
theory has kept the SU(4) symmetry of the model. The 
present results are in conflict with the prediction of a 
spin liquid state with a finite energy gap by the quan- 
tum Monte Carlo Green function calculations on a small 
square lattice for the same model [0. However, the very 
small magnetization clearly indicates that the quantum 
fluctuations in the spin-orbital coupled model are much 
stronger than those in the pure spin model. The numer- 
ical data shows that the energy gap (even if it exists) 
is very close to zero. In this sense, both analytical and 
numerical results seem to be consistent with each other. 

Several remarks are in order, i) The reduced model 
Hamiltonian has displayed a unique form of attractive 
pairing instability among the nearest neighbor hard-core 
bosons. Introducing a short-ranged VB pairing order pa- 
rameter does not break the global SU(4) symmetry of the 
model, ii) The mean field theory in ID has given rise to 



the correct ground state of a quantum disordered spin- 
orbital liquid state, which was consistent with numerical 
results ||l^. iii) On a 2D square lattice, the long-range 
spin-orbital ferromagnetic ordered state is obtained due 
to the Bose-Einstein condensation. We can also expect 
that such a long-ranged ordered state should be robust to 
the perturbations added even away from the SU(4) sym- 
metric point slightly, iv) Compared with the numerical 
method, the SU(4) Schwinger boson mean field approach 
for the symmetrically coupled spin-orbital systems has 
provided the insight on the exact nature of the ground 
state and its first excitations. 

In conclusion, we derive a reduced effective spin-orbital 
Hamiltonian in a symmetric spin-orbital systems by in- 
troducing the SU(4) Schwinger boson representation. Af- 
ter introducing a symmetric pairing VB order parame- 
ter, a mean field theory has been developed, leading to a 
quantum disordered liquid state with a finite gap in the 
collective excitations in ID, and on a 2D square lattice 
a novel long-range FM ordered spin-orbital crystal state 
where both the spin and orbital degrees of freedom form 
AFM. The new ordering properties due to the interplay 
between the spin and orbital degrees of freedom can be 
detected in the future experiments on transitional metal 
oxides. 
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